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It is shown that dynamical screening of the Coulomb interaction in LaCuO leads to low-energy 
electronic collective excitations in a small region around the c-axis, strongly mixing with certain 
interlayer phonons polarized along this axis. The manifestation of such a phonon-plasmon scenario 
in layered systems based on a nonadiabatic charge response is quantitatively supported by a realistic 
calculation of the frequency and wavevector dependent irreducible polarization part of the density 
response function (DRF). The latter is used within linear response theory to calculate the coupled 
mode dispersion in the main symmetry directions of the Brillouin zone (BZ) and the charge density 
redistributions excited by certain strongly coupling phonon-like and plasmon-like modes. Moreover, 
the corresponding mode induced orbital averaged changes of the selfconsistent potential felt by the 
electrons are assessed. Our analysis should be representative for the optimally to overdoped state 
of the cuprates where experimental evidence of a coherent three-dimensional Fermi surface and a 
coherent c-axis charge transport is given. It is demonstrated that modes from the outside of a 
small region around the c-axis can reliably be calculated within the adiabatic limit. Only a minor 
nonadiabatic correction is found for these modes in form of stiffening of the high-frequency oxygen- 
bond-stretching mode at the X-point which is attributed to dynamical reduced nesting. On the 
other hand, modes inside the nonadiabatic sector of the BZ have to be determined nonadiabatically 
owing to the poor dynamical screening of the long-ranged Coulomb interaction around the c-axis 
by the slow charge dynamics. In particular, the relevance of the strongly coupling phonon-like apex 
oxygen Z-point breathing mode at about 40 meV is emphasized whose mode energy decreases with 
less doping. Comparing the calculations with experimental neutron scattering and infrared mea- 
surements provides strong evidence for low-lying plasmons around the c-axis coupling strongly with 
corresponding optical phonons. The existence of these coupled modes evidenced by the calculations 
yields additional virtual exchange bosons for pairing. Finally, the strong nonlocal, nonadiabatic 
polar electron-phonon coupling around the c-axis found in the computations means that besides the 
repulsive short-ranged part of the Coulomb interaction responsible for strong correlation effects also 
the long-ranged part plays an important role for the physics in the cuprates. 
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I. INTRODUCTION 

In layered, strongly anisotropic materials like the cop- 
per oxide high-temperature superconductors (HTSC's) 
the screening of the Coulomb interaction between the 
layers is imperfect and its dynamic nature becomes im- 
portant. Perpendicular to the CuO planes low- lying elec- 
tronic collective modes (interlayer plasmons) can be ex- 
pected. In connection with superconductivity the inter- 
play between the attractive interaction, mediated e.g. 
by phonons, and the dynamically screened Coulomb 
interaction is more complex than in the conventional 
theory. Accordingly, the existence and possible rele- 
vance of energetically low lying plasmons, e.g. accous- 
tic plasmons, for superconductivity in the cuprates has 
been discussed before in the literature, see e.g. Refs. 

5.6. A direct observation of low lying plasmons 
in the HTSC's, e.g. by high-resolution inelastic elec- 
tron scattering or some other technique is still missing. 
As shown in this work there is only a small q-space re- 
gion around the c-axis with a nonadiabatic charge re- 
sponse where low lying phonon-plasmon modes can ex- 
ist and thus a very high q-space resolution would be 



needed in an experiment. However, there is indirect ex- 
perimental evidence from different probes, like the use 
of reflectance and ellipsometric measurements 7 -, external 
losses in photoemission 8 , infrared reflectivity^ or reso- 
nant inelastic ir-ray scattering^. 

A proper account of dynamical screening of the 
Coulomb interaction is not only substantial for super- 
conductivity, but also for a correct description of normal 
state properties of the HTSC's such as phonon dynam- 
ics, see our simplified model approach in Refs. [Tllfl3 and 
references therein. These calculations indicate that in a 
small sector around the c-axis of LaCuO phonon-plasmon 
mixing becomes likely if the interlayer coupling is suffi- 
ciently weak. The reason for such a mixing behaviour is a 
sufficiently slow electron dynamics perpendicular to the 
CuO planes. However, the important question remains 
if this is really true in LaCuO or other cuprates. This 
topic and its consequences is investigated in the present 
work. Electron dynamics and phonon dynamics then 
have about the same time scale and as a consequence 
nonadiabatic behaviour results. 

Due to the incomplete dynamical screening of the 
Coulomb interaction a strong nonlocal, polar electron- 
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phonon coupling results along the c-axis together with 
optical conductivity, also in the well-doped metallic state. 
Experimentally this can be seen for example by the in- 
frared measurements for LaCu 13 ' 14 i 15 ' 16 . Direct evi- 
dence for the presence and importance of this type of un- 
conventional electron-phonon coupling around the c-axis 
where the electrons in the cuprates are strongly coupled 
to c-axis polar phonons also comes from a recent study 
of the self-energy of the nodal quasiparticles in Bi2201 
(Ref. [I?]). For a discussion of phonon-plasmon mixing 
in the cuprates in the context of many-body polaronic 
effects in the phonon spectrum, we refer to Ref. 

As far as underdoped and even optimally doped 
LaCuO is concerned, the optical c-axis spectra display 
the features typical for an ionic insulatori 3 ^^^ and 
not those of a metal dealt within adiabatic approxima- 
tion. The spectra are dominated by optical phonons and 
are almost unchanged from that of the insulating parent 
compound upon doping. These experimental facts can- 
not be explained by the first principles calculations of 
the phonon dynamics of the HTSC's published so far in 
the literatur o 19 ' 20 i 21 ' 22 because these computations have 
been performed within (static) density functional the- 
ory (DFT) mostly in local density approximation (LDA) . 
Such calculations are based on the adiabatic approxi- 
mation and always yield a static metallic screening be- 
haviour consistent with closed LO-TO splittings, at the 
T-point, i.e. in particular closed ^^-splittings in LaCuO. 
This is clearly in contrast to the optical activity in the 
infrared experiments. In static DFT calculations the 
transverse effective charges vanish in the metallic state 
and consequently the induced dipole moments defining 
the oscillator strenghts in the dielectric function (matrix) 
vanish, too. So, there will be no optical activity by the 
phonons. 

In the first principles calculation o 19 ' 20 ^ 21 ' 22 also no ex- 
plicit results for the phonon dispersion curves along the 
c-axis or a small region around this axis have been pre- 
sented. However, this has been accomplished within a 
simplified model in Refs. [lllfl2l and will quantivatively 
be computed in this work. In Ref. [13 only a short re- 
mark is given for the case of LaCuO stating that certain 
experimental results along the c-axis are not supported 
by the calculation. Moreover, when judging the findings 
in Refs. I22II23I of a small phonon contribution to the 
photoemission kink detected by ARPES 2 -i2£ one should 
keep in mind that these very interesting calculations also 
have been performed within static DFT, neglecting e.g. 
the strong nonadiabatic dynamically screened polar cou- 
pling around the c-axis and also correlation effects be- 
yond LDA. On the other hand, the importance of the 
electron-phonon interaction in context with the kink in 
the dispersion of the nodal quasiparticles is emphasized 
in Ref. [26| on the basis of the two-dimensional three- 
band Hubbard model with electron-phonon interaction 
included and also in Ref. H3. 

Because of the strong anisotropy of the cuprates and 
consequently an expected very weak dispersion of the 



electronic bandstructure along the c-axis, which, how- 
ever, is systematically overestimated by DFT-LDA calcu- 
lations, electron dynamics and c-axis phonon dynamics is 
coupled dynamically and a nonadiabatic treatment with 
a dynamical screened Coulomb interaction is essential. 
On the other hand, first principles calculation o 19 i 20 ' 21 i 22 
and also our microscopic calculations, see e.g. Refs. 
[l2ll28ll29l . have shown that the static DFT is in general a 
sufficient approximation for phonon modes propagating 
in the CuO-plane of the HTSC's, but probably not so 
around the c-axis. 

For modes from this q-space region a coupled phonon- 
plasmon scenario has been suggested for LaCuO within 
a simplified microscopic model approach in the frame- 
work of linear response theor y 11 ! 30 assuming a sufficiently 
weak interlayer coupling. In this treatment an eleven- 
band model (11BM) that represents the two-dimensional 
electronic structure of the CuO plane in terms of the 
Cu3d and O xy 2p orbitals is generalized to the third di- 
mension by introducing a suitable interlayer coupling 
in parametrized form. By varying this coupling from 
the outside it is possible to pass from the strictly two- 
dimensional nonadiabatic case of the charge response 
(charge confinement in the CuO plane) to a moderate 
anisotropic adiabatic situation, typical for DFT-LDA 
calculations. A comparison with experimental results 
from inelastic neutron scattering (INS) and infrared spec- 
troscopy also has been given. From the interpretation of 
the results a phonon-plasmon mixing seems to be likely. 
Adopting within such a model approach a suitable cou- 
pling we were able to explain at least qualitatively an 
appearent inconsistency between the current interpreta- 
tion of the INS results for the Ai ~ (0, 0,1) branches 
in doped metallic samples of LaCuO (Refs. [31. 32) . The 
branches look as expected in adiabatic approximation for 
a less anisotropic metal, i.e. featuring closed A% u split- 
tings. On the other hand, the infrared data display as 
already mentioned significant ^u-sphttings typical for 
an ionic insulator, even for well doped metallic probes. 
Moreover, in Ref. [ll] arguments have been presented in 
the framework of phonon-plasmon mixing to understand 
at least qualitatively the large softening and the massive 
line broadening of the apical oxygen breathing mode Of, 
at the Z-point of the BZ. Altogether, the model calcula- 
tions point to phonon-plasmon coupling in a small nona- 
diabatic region around the c-axis where the Coulomb in- 
teraction has to be screened dynamically. A direct ex- 
perimental search for such a scenario is inconclusive so 
fai— — . 

From our calculations in this work we can conclude 
that much better q-space resolution transverse to the c- 
axis would be needed in the experiment to resolve the 
coupled mode dispersion in the small nonadiabatic sector. 

One important topic of this paper is to scrutinize and 
quantify the discussion performed within the simplified 
model for LaCuO by starting with a realistic tight bind- 
ing approximation (TBA) of the first principles band- 
structure (31-band-model (31BM), Ref. l3~il ) which is as 
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a typical LDA-like bandstructure not anisotropic enough 
and adapt the latter to the real anisotropy of the ma- 
terial. With such a realistic electronic bandstructure 
the wavevector and frequency-dependent irreducible po- 
larization part of the DRF is calculated. The latter is 
used in linear response theory to obtain the dynami- 
cally screened Coulomb interaction and ultimately the 
full nonadiabatic coupled mode dispersion in the main 
symmetry directions of the BZ. Moreover, charge density 
redistributions Sp of certain strongly coupling phonon- 
like and plasmon-like modes of the nonadiabatic region 
and the corresponding mode induced changes of the self- 
consistent potential SV felt by the electrons are inves- 
tigated and compared with the adiabatic limit. From 
our calculations a phonon-plasmon mixing and a strong 
nonlocal, nonadiabatic polar electron-phonon interac- 
tion is predicted for LaCuO. This conclusion should be 
representative for the optimally to overdoped state of 
the cuprates where experimental evidence of a coherent 
three-dimensional Fermi surface (FS) 3 ^ and coherent c- 
axis charge transport is give n 36 ' 37 ! 38 . 

The article is organized as follows. In Sec. II some el- 
ements of the theory necessary to understand the calcu- 
lated results are reviewed. Section III presents the com- 
putations. It gives a discussion of the electronic struc- 
ture and the nonadiabatic and adiabatic mode dispersion. 
Furthermore, Sp and SV are investigated for certain rel- 
evant generic modes in the cuprates. Finally, a summary 
of the paper is presented in Sec. IV and the conclusions 
are drawn. 



II. ELEMENTS OF THE THEORY AND 
MODELING 

In the following a survey of the theory and modeling 
is presented. A detailed description can be found in Ref. 
l39l and in particular in Ref. 40 where the calculation of 
the coupling parameters of the theory is presented. 

The local part of the electronic charge response and 
the EPI is approximated in the spirit of the quasi-ion 
approach 2 *^ by an ab initio rigid ion model (RIM) 
taking into account covalent ion softening in terms of 
(static) effective ionic charges calculated from a tight- 
binding analysis. The tight-binding analysis supplies 
these charges as extracted from the orbital occupation 
numbers Q M of the p (tight-binding) orbital in question: 

^ = §£|cW(k)l 2 - (1) 

nk 

Cfm (k) stands for the /Lt-component of the eigenvector of 
band n at the wavevector k in the first BZ; the summa- 
tion in (TJJ) runs over all occupied states and N gives the 
number of the elementary cells in the (periodic) crystal. 

In addition, scaling of the short-ranged part of certain 
pair potentials between the ions is performed to simulate 
further covalence effects in the calculation in such a way 



that the energy-minimized structure is as close as pos- 
sible to the experimental one^ 2 -. Structure optimization 
and energy minimization is very important for a reliable 
calculation of the phonon dynamics through the dynam- 
ical matrix. Taking just the experimental structure data 
as is done in many cases in the literature may lead to 
uncontrolled errors in the phonon calculations. 

The RIM with the corrections just mentioned then 
serves as an unbiased reference system for the descrip- 
tion of the HTSC's and can be considered as a first ap- 
proximation for the insulating state of these compounds. 
Starting with such an unprejudiced rigid reference sys- 
tem non-rigid electronic polarization processes are intro- 
duced in form of more or less localized electronic charge- 
fluctuations (CF's) at the outer shells of the ions. Es- 
pecially in the metallic state of the HTSC's the latter 
dominate the nonlocal contribution of the electronic den- 
sity response and the EPI and are particularly impor- 
tant in the CuO planes. In addition, anisotropic dipole- 
fluctuations (DF's) are admitted in our approac h 40 ' 43 , 
which prove to be specifically of interest for the ions in 
the ionic layers mediating the dielectric coupling and for 
the polar modes. Thus, the basic variable of our model 
is the ionic density which is given in the perturbed state 
by 

p a (r,Qx,p a ) = p ( Ur) + J2QxPx F (r) + Pa-rpV(r). (2) 

A 

p a is the density of the unperturbed ion, as used in the 
RIM, localized at the sublattice a of the crystal and mov- 
ing rigidly with the latter under displacement. The Q\ 
and p^ F describe the amplitudes and the form-factors of 
the CF's and the last term in Eq. ([2|) represents the dipo- 
lar deformation of an ion a with amplitude (dipole mo- 
ment) p a and a radial density distribution p® . f denotes 
the unit vector in the direction of r. The p^ F are approx- 
imated by a spherical average of the orbital densities of 
the ionic shells calculated in LDA taking self-interaction 
effects (SIC) into account. The dipole density p® is ob- 
tained from a modified Sternheimer method in the frame- 
work of LDA-SIC 40 . All SIC-calculations are performed 
for the average spherical shell in the orbital-averaged 
form according to Perdew and Zunger— . For the correla- 
tion part of the energy per electron e the parametrization 
given in Ref. [44| has been used. 

The total energy of the crystal is obtained by assum- 
ing that the density can be approximated by a superpo- 
sition of overlapping densities p a . The p a in Eq. ^ 
are also calculated within LDA-SIC taking environment 
effects, via a Watson sphere potential and the calculated 
static effective charges of the ions into account. The 
Watson sphere method is only used for the oxygen ions 
and the depth of the Watson sphere potential is set as 
the Madelung potential at the corresponding site. Such 
an approximation holds well in the HTSC As a 

general rule, partial covalence reduces the amplitude of 
the static effective charges in mixed ionic-covalent com- 
pounds like the HTSC's, because the charge transfer from 
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the cations to the anions is not complete as in the entirely 
ionic case. Finally, applying the pair-potential approxi- 
mation we get for the total energy: 

e(r,o=y, e *(o + \ E ^(r^-rs.c). 

a,a (a.a)^(b,/3) 

(3) 

The energy E depends on both the configuration of the 
ions {R} and the electronic (charge) degrees of freedom 
(EDF) {C} of the charge density, i.e. {Q\} and {p a } in 
Eq. ([2]). E* are the energies of the single ions, a, b 
denote the elementary cells and a, (3 the corresponding 
sublattices. The second term in Eq. ^ is the interaction 
energy of the system, expressed in terms of anisotropic 
pair-interactions $> a p- Both E* and & a p in general de- 
pend upon £ via p a in Eq. @. 

The pair potentials in Eq. ([3]) can be seperated into 
long-ranged Coulomb contributions and short-ranged 
terms as follows: 



with 



$a/3(R, C) 



Z a Zp 
R 



- (ZaPp + ZpPa) 



i? 3 



Pa ' P/3 



R 3 



3 (p a .R)(R. P/3 ) + ^ (R;a 



<E> a/3 (R, C) =K a Up(R, C) + KpU a (R, C) 
+ W Q/3 (R,C) + G^(R,C). 



(4) 



(5) 



The first term in Eq. ([4]) describes the long-ranged ion- 
ion, the second the dipole-ion and the third and fourth 
term the dipole-dipole interaction. Z a and Zp are the 
variable charges of the ions in case the CF's are excited. 
The latter reduce to the ionic charges for rigid ions. K a 
and Kp are the charges of the ion cores. The remain- 
ing term in Eq. (|H) given in ([5]) represents the short- 
ranged interactions. These short-ranged contributions to 
the pair potentials are expressed by the following inte- 
grals: 



U a (R,() =- J d 3 rp a (r, 



C)x 



Ir-R| 



1 

R 



r R 



(6) 



W af3 (R,0 = J d 3 r J dV[p Q (r,C)Mr',Ox 
1 1 (r + r') • R 



|r - r' - R| R 



R :i 



, (?) 



G a0 (TL,O= J d 3 r[p a p(r,()e(p a p(r,()) 



Pa(r,C)e(Po(r,C)) 
^(r-R,C)e(p/3(r-R,C))], 



(8) 



p a p(r, C) = Pa(r, C) + Pp(r - R, C)- 



(9) 



K a Up(R., £) yields the short-ranged contribution of the 
interaction between the core a and the density pp accord- 
ing to Eq. ([2]) . W a p (R, £) represents the short-ranged 
Coulomb contribution of the interaction of the density 
p a with the density pp. G a p(R,£) is the sum of the 
kinetic one-particle- and the exchange-correlation (XC) 
contribution of the interaction between the two ions^S. 
The short-ranged part of the potentials and the various 
coupling coefficients are calculated numerically for a set 
of disctances R between the ions. The corresponding re- 
sults are than described by an analytical function of the 
form: 



f(R) = ±cxp (a + 0R- 



(10) 



a, [3 and 7 in Eq. (|10[) are fit parameters. From the 
adiabatic condition 



dE(R,Q 







(11) 



an expression for the atomic force constants, and accord- 
ingly the dynamical matrix in harmonic approximation 
can be derived: 



^(q) 



RIM 



(12) 

The first term on the right hand side denotes the contri- 
bution from the RIM. M a , Mp are the masses of the ions 
and q is a wave vector from the first BZ. 

The quantities B(q) and C(q) in Eq. (TT2"]) represent 
the Fourier transforms of the electronic coupling coeffi- 
cients as calculated from the energy in Eq. ([3]), or the 
pair potentials in Eqs. (J4j)— (J9j> , respectivly. 



Bab 
K P 



Cab 
u?ur' 



d 2 E(R,Q 

d(*dR b p ■ 

d 2 E(R,Q 



(13) 
(14) 



k denotes the EDF (CF and DF in the present model, 
see Eq. ©) in an elementary cell. The B coefficients 
describe the coupling between the EDF and the dis- 
placed ions (bare electron-phonon coupling), and the co- 
efficients C determine the interaction between the EDF. 
The phonon frequencies Wo-(q) and the corresponding 
eigenvectors e a (qa) of the modes (qcr) are obtained from 
the secular equation for the dynamical matrix in Eq. 
CD, i.e. 



E 

0j 



^(q)ef(q)=c 2 (q) e ?(q). 



(15) 
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The Eqs. ([T2|) - ([T5|) are generally valid and, in particular, 
are independent of the specific model for the decompo- 
sition of the perturbed density in Eq. ((2|) and the pair 
approximation Eq. ([3]) for the energy. The lenghty de- 
tails of the calculation of the coupling coefficients B and 
C cannot be reviewed in this paper. They are presented 
in Ref. [40t In this context we remark that the coupling 
matrix C KK '(q) of the EDF-EDF interaction, whose in- 
verse appears in Eq. (| 1 2() for the dynamical matrix, can 
be written in matrix notation as 



C = IT 



V. 



(16) 



IT" 1 is the inverse of the irreducible polarization part 
of the density response function (matrix) and contains 
the kinetic part to the interaction C while V embodies 
the Hartree and exchange-correlation contribution. C _1 
needed for the dynamical matrix and the EPI is closely 
related to the (linear) density response function (matrix) 
and to the inverse dielectric function (matrix) e -1 , re- 
spectively. 

Only very few attempts have been made to calculate 
the phonon dispersion and the EPI of the HTSC's using 
the linear response method in form of density functional 
perturbation theory (DFPT) within LDAi^ 2 ^. These 
calculations correspond to calculating II and V in DFT- 
LDA and for the metallic state only. On the other hand, 
in our microscopic modeling DFT-LDA-SIC calculations 
are performed for the various densities in Eq. ^ in order 
to obtain the coupling coefficients B and V. Including 
SIC is particularly important for localized orbitals like 
Cu3d in the HTSC's. SIC as a correction for a single 
particle term is not a correlation effect, which per def- 
inition cannot be described in a single particle theory, 
but SIC is important for contracting in particular the 
localized Cu3d orbitals. Our theoretical results for the 
phonon dispersio n 12 i 28 ' 43 , which compare well with the 
experiments, demonstrate that the approximative calcu- 
lation of the coupling coefficients in our approach is suf- 
ficient, even for the localized Cu3d states. Written in 
matrix notation we get for the density response matrix 
the relation 



c- 1 = 11(1 



Vtt)- 1 = lie- 1 , 



e = 1 + VH. (17) 



The CF-CF submatrix of the matrix II can approxima- 
tively be calculated for the metallic (but not for the un- 
doped and underdoped) state of the HTSC's from a TBA 
of a single particle electronic bandstructure. In this case 
the electronic polarizability n in tight-binding represen- 
tation reads: 

n -'(^- = °) = -]v 2J k ^ (k + q) -^( k) x 

x [C*„(k)C K „,(k + q)] [C:,„(k)CW(k + q)]* . 

(18) 

/, E and C in Eq. (fT8|) are the occupation numbers, the 
single-particle energies and the expansion coefficientes of 
the Bloch-functions in terms of tight-binding functions. 



The self-consistent change of an EDF at an ion induced 
by a phonon mode (qcr) with frequency w CT (q) and eigen- 
vector e°(qcr) can be derived in the form 



-]rX™(q)u Q ( qf 7) 



with the displacement of the ions 
<(qcr) = 



(19) 



1/2 



2M Q w CT (q) 



e a (qcr)e lqRa 



(20) 



The self-consistent response per unit displacement of the 
EDF in Eq. (|T9^) is calculated in linear response theory 
as: 



X(q) = n(q)e- 1 (q)B(q) = C- 1 (q)B(q). 



(21) 



A measure of the strength of the EPI for a certain phonon 
mode (qcr) is provided by the change of the self-consistent 
potential in the crystal felt by an electron at some space 
point r, i.e. 6V e ff{r, qcr). Averaging this quantity with 
the corresponding density form factor p K (r — R*) at the 
EDF located at R* we obtain 



SV^qa) 



dVp K (r-B.*)6V g{r,qa) 



(22) 



as a parameter for the strength of the EPI in the mode 
(qcr) mediated by the EDF considered. For an expression 
of 6V*(qa) in terms of the coupling coefficients in Eqs. 
(fI5]) and ([H]). see Ref. QJ. 

The generalization for the quantity n in Eqs. (p~6|) and 
(jTTJ) needed for the kinetic part of the charge response 
in the nonadiabatic regime, where dynamical screening 
effects must be considered, can be achieved by adding 
(hu + in) to the differences of the single-particle ener- 
gies in the denominator of the expression for n in Eq. 
(jT5J). Other possible nonadiabatic contributions to C re- 
lated to dynamical exchange-correlation effects and the 
phonons themselves are beyond the scope of the present 
model. Using Eq. (|17p for the dielectric matrix, e, and 
the frequency-dependent version of the irreducible polar- 
ization part, n, according to Eq. (TT8"]) . the free-plasmon 
dispersion is obtained from the condition, 



det[e reK /(q,w)] = 0. 



(23) 



The coupled-mode frequencies of the phonons and the 
plasmons must be determined self-consistently from the 
secular equation (|15p for the dynamical matrix which now 
contains the frequency u implicitly via n in the response 
function C -1 . Analogously, the dependence on the fre- 
quency is transferred to the quantity X in Eq. (|21[) and 
thus to 6( K and SV K in Eqs. (| 19[) and ([22]) , respectively. 
Such a nonadiabatic approach is necessary for a descrip- 
tion of the interlayer phonons and the charge-response 
within a small region around the c-axis. 
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FIG. 1: Electronic bandstructure E n (\s) of LaCuO in the 31 
band model (31BM) obtained from an accurate tight-binding 
representation of the first-principles linearized-augmented - 
plane- wave bandstructure^. The Fermi energy Ef corre- 
sponding to the undoped material is taken as the zero of the 
energy (a). Density of states Z(E) of the 31BM (b). 



The time-consuming numerical calculations were car- 
ried out on the computers of the Morfeus GRID at the 
Westfalische Wilhelms-Universitat Miinster, with the use 
of Condor—. 



III. RESULTS AND DISCUSSION 

A. Model of the single particle content of the 
irreducible polarization part 

As a first approximation the electronic bandstructure 
(BS) E n (k) of LaCuO and the expansion coefficients 
C ren (k), needed as input for the single particle content 
of the irreducible polarization part H KK > in Eq. (|18|) 
are taken from an accurate tight-binding representation 
of the first-principles linearized-augmented-plane-wave 
bandstructure (LAPW) as obtained within the frame- 
work of DFT-LDA21 

This analysis leads to a 31BM including the La5d, 
Cu3d, 4s, 4p, and 02p states. The result for the re- 
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FIG. 2: Calculated phonon dispersion for LaCuO in the adi- 
abatic approximation for the c-axis polarized Ai modes (~ 
(0,0,1)) based on the bandstructure within the 31BM as input 
for the proper polarization part II KK / (a), (b)-(d) calculated 
Ai modes as in (a) but using the modified 31BM (M31BM) 
as discussed in the text as input for TL KK i. The panels (b)- 
(d) represent the phonon dispersion for three different doping 
levels x = 0.156 (E F = -0.08), x = 0.224 (E F = -0.10) and 
x — 0.297 (Ef = —0.135), by applying a rigid-band approx- 
imation of the M31BM to represent the optimally doped to 
overdoped state of LaCuO. The exper imental results (open 



expe 

squares □) are taken from Refs. [31 32. The full dot (•) de- 
notes the Of mode. Ef is in units of eV. 



calculated tight-binding bandstructure is shown in Fig. 
HJa) and the corresponding density of states (DOS) 



(24) 



is displayed in Fig. QJb). The van Hove peak in the DOS 
is broadened by the dispersion of the bandstructure in 
^-direction. 

The corresponding calculated phonon dispersion of 
LaCuO in the adiabatic approximation for the Ai modes 
polarized along the c-axis being most sensitive with re- 
spect to the interlayer coupling is presented in Fig. H^a). 
The open squares indicate the present " interpretation" in 
the literature of the experimental result a 31 i 32 . This inter- 
pretation will be discussed in detail in Sec. Ill B. The 
characteristic experimental features of the dispersion are 
the step-like structure of the second highest branch and 
most significant the third highest branch with the steep 
dispersion from the T-point towards the Z-point. Both 
features are not well reflected in the calculation based on 
the typical DFT-LDA-like bandstructure of the 31BM 
which underlies the computation of the static (u> = 0) 
proper polarization part Yl KK > in Eq. (JTHj). As we see in 
a moment, the reason for this is an overestimation of the 
residual k z dispersion of the bands, i.e. the DFT-LDA- 
like bandstructure is not anisotropic enough. It should be 
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remarked that an earlier calculation within the 31BM2£ 
neglecting dipole fluctuations as additional electronic po- 
larization processes besides the charge fluctuations still 
further enhances the deviations between the calculated 
and the measured result for the Ai-modes. Anisotropic 
DF's being particular important along the c-axis and for 
ions in the ionic layers^ have been taken into account 
in the calculations shown in Fig. [DJa) and in all other 
computations of the phonon dispersion presented in this 
paper. 

In order to improve the LDA-like electronic BS of the 
3 IBM which overestimates the c-axis coupling we have 
investigated the effect on the Ai-modes of a modification 
of the tight-binding parameters of the 31BM. It turns out 
that a reduction of the first neighbour O^-La parame- 
ters by 1/6 and of the first neighbour La-La parameters 
by 1/3 leads to a much better result for the calculated 
Ai-modes in adiabatic approximation, see Fig. 2(b-d). 
The latter results have been obtained for three different 
doping levels, where the effect of alloying was treated 
in rigid-band approximation by lowering the Fermi level 
appropriately to accommodate x holes per primitve cell. 

The modified electronic BS of LaCuO and the corre- 
sponding DOS is shown in Fig. [3ja) and Fig. 02(b), 
respectively. The calculated phonon dispersion is now 
in good agreement with the experiment. In particular 
the characteristic features, i.e. the Ai-branch with the 
steep dispersion and the step-like behaviour are well de- 
scribed. Moreover, we find as a consequence of the en- 
hanced anisotropy a rearrangement of the three Z-point 
modes with the highest frequencies. While in the 31BM 
the anomalous Of mode is the second highest mode in 
the more anisotropic modified 31BM (M31BM) with a 
weaker coupling along the c-axis Of is the lowest of the 
three modes and the end point of the steep branch. 

From a comparison of the BS of the 31BM in Fig. [TJa) 
with the BS of the M31BM in Fig 3(a) we extract that 
in the latter case the saddle point region around A/2 
is more extended and the width of the BS is reduced 
by about 0.5 eV in the energy range shown in the fig- 
ure. Important for the charge response along the c-axis 
is the decrease of the /c z -dispersion of the bands along 
the A direction, which obviously is consistent with the 
experimental phonon dispersion of the Ai modes being 
most sensitive to the charge dynamics along the c-axis. 
Comparing the calculated results for the DOS in both 
models, see Fig. [Tfb) and Fig. [3jb) , respectively, we find 
that the broadening of the van Hove peak is decreased in 
the M31BM by about 50 % because of the reduced k z - 
dispersion. Simultaneously, the DOS Z(E) around the 
peak is enhanced by the increased anisotropy. 

Another important effect of the weakening of the in- 
terlaycr coupling in our modeling as compared to a typi- 
cal DFT-LDA bandstructure, leading to a significant en- 
hancement of the DOS Z(Ep) at the Fermi energy Ep, is 
the amplification of the proper polarization part Tl KK i (q) 
on his part. This can be seen in the long wavelength limit 
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FIG. 3: Electronic bandstructure E n (k) of LaCuO in the 
modified 31 band model (M31BM) taking into account the 
enhanced anisotropy of the real material as compared with 
the 3IBM typical for LDA (a). Density of states Z(E) in the 
M31BM (b). The vertical lines indicate the different doping 
levels from optimal to overdoping as in Fig. [2] 



where the sum rule 



lim^n KK ,(q) = .Z(£ F ) 



(25) 



rigorously holds for a metal^. Altogether, the re- 
sults demonstrate the importance of a correct interlayer- 
coupling also for the electronic properties in the CuO 
layer. 

The vertical lines in Fig. 03b) indicate three doping 
levels x in rigid-band approximation. The specific dop- 
ings have been selected to model the optimally doped 
state of LaCuO with the Fermi energy Ep at the van 
Hove peak (model OP: E F = -0.08 eV, x = 0.156) 
and two overdoped states (model ODI: Ep — —0.10 eV, 
x = 0.224; model OD2: Ep = -0.135 eV, x = 0.297). 

The doping levels chosen in this way allow for a 
comparison of measured Fermi surfaces (FS) at these 
levels^ii^ w ith the calculated FS within the M31BM, see 
Fig. |U From the comparison we conclude that the overall 
features of the BS calculation in the M31BM are in good 
agreement with the results obtained by angle-resolved 
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FIG. 4: Comparison of the doping evolution of the measured 
Fermi surface of La,2- x Sr x CuCU from Ref. [4^(4^ for optimal 
doping (x = 0.15) to overdoping (x — 0.22; x — 0.30) with 
the corresponding calculated Fermi surface within the rigid- 
band approximation of the M31BM, respectively. The dots 
indicate the experimental results. 



photoemission spectroscopy (ARPES) over the doping 
range considered. Such an agreement obtained within 
the M31BM and the consistence with the Ai phonons 
speaks in favour of this model with renormalized c-axis 
hopping parameters leading to an enhanced anisotropy. 
Moreover, it demonstates the importance of an experi- 
mental verification of a computed electronic bandstruc- 
ture. Clearly visible in the experiments and the calcula- 
tions is the change of the FS topology and of the related 
nesting structures upon doping when passing from the 
model OP to model OD1 or OD2, respectively. 

Strong nesting of the FS can bring about either a 
charge density wave or a spin density wave, or possibly 
both. The tendency toward Fermi surface driven insta- 
bilities frequently is anounced by maxima of the so called 
(noninteracting) suszeptibility 



n (q) 



- V 

N ^ 



/ n ,(k + q)-/ n (k) 
£ n ,(k + q)-£„(k)' 



(26) 



which can be obtained form Eq. (fl8|) by equating all the 
expansion coefficients C Kn (k) to one. It is obvious from 
this form that possible FS-nesting at certain q vectors is 
reflected as a maximum of this function. 

From inspection of the FS in Fig. 3] we can expect 
that approximate nesting vectors around the X-point of 
the BZ will lead to corresponding maxima in llo (q) . For 
the optimally case (model OP) we find a maximum at 
the X-point and in the overdoped state (model OD1 
and OD2) the maximum is progressively shifted toward 
the r point. The position of the maximum for model 
OP is significant because it coincides with the ordering 
wave vector in the antiferromagnetic state. For the over- 
doped state, the ordering wave vector (magnetic scat- 
tering peak) is predicted from our calculations for llo 
to move inward from the X-point. Quite general, in the 
doped cuprates the spin fluctuations (SF's) are antiferro- 
magnetic in origin and arise from both, nesting properties 
of the FS (spin-density-wave type SF's) and from cor- 
relation driven nearest-neighbour antiferromagnetic su- 
perexchange because of the proximity of a long-range 
ordered antiferromagnetic state in the undoped parent 



compounds. 

In Fig. [5]we display for the models OP, OD1 and OD2 
FS maps for different k z values (k z = 0, k z — 2zl) as a 
function of k = (k x , k y ). This allows for example to map 
the k x -, fcy-dependence of the c-axis dispersion. The cal- 
culations demonstrate that the effect of the k z dispersion 
vanishes along the nodal direction and that it increases 
as one moves towards the antinodal region. The three- 
dimensionality associated with the £; z -dispersion, usually 
neglected in discussing the physics of the cuprates, has 
been shown to play a key role in shaping the ARPES 
spectr a 36 ! 50 , because the residual A: z -dispersion of bands 
in a quasi-two-dimensional material will induce an irre- 
ducible linewidth in ARPES peaks. On the other hand 
detecting such k z related linewidth in the ARPES spec- 
tra as found for LaCu02S establishes the existence of 
coherent c-axis transport. 

Finally to get also a more global impression of the mag- 
nitude of the enhanced anisotropy in the M31BM as com- 
pared with a typical DFT-LDA based first principles re- 
sult we compare some FS parameters being important for 
transport properties like the Drude plasma energy ten- 
sor and the Fermi velocity tensor. The Drude tensor is 
defined as 



HQ 



8tt x -, 

NV z j^ 

^ kn 



5{E n (k)~E F ) 



and the Fermi velocity tensor ist given by 



with 



Vkr, 



1 dE n (k) 
h <9k ' 



(27) 



(28) 



(29) 



The outcome from our model OP for the optimally doped 
state is compared in Table U with the corresponding 
results calculated for optimally doped LaCuO within 
LAP W 51 ' 52 . From the table we extract an enhance- 
ment of the anisotropy ratio A 



- - n 



/Op 
i a facto: 

M31BM as compared to the LAPW calculation 



and of 



A VF = {v 2 F xx ) 1 ^ 2 1 {v 2 F zz ) 1 ^ 2 by about a factor 5 in the 



B. Nonadiabation phonon calculations and 
phonon-plasmon coupling 

We prelude this section with a discussion of the anoma- 
lous behaviour of the apex-oxygen bond-stretching mode 
in LaCuO at the Z-point of the BZ (Of) polarized per- 
pendicular to the CuO plane, see Fig. \§\ From the ex- 
perimental side, it took quite some time until this mode 
could be assessed^ 2 - because of a massive line-broadening 
of about 4 THz in the metallic phase at optimal doping 
and its very large softening of about 5.5 THz when pass- 
ing from the insulating to the metallic stat o 31 i 32 . The 



(a) (b) (c) 

FIG. 5: Fermi surface maps for different k z values (k z = z = and k z = z = — ) as a function of k|| = (kx,k y ) for different 
doping levels in rigid-band approximation of the M31BM. Model OP (x = 0.156) (a), model OD1 (x = 0.224) (b) and model 
OD2 (x = 0.297) (c). 













An 




LAPW 


701.21 


132.99 


2.20 


0.41 


5.27 


5.37 


OP 


648.60 


25.25 


2.99 


0.11 


25.69 


27.97 



Table I: Comparison of the Fermi surface parameters (Drude 
plasma energy tensor, Fermi velocity) and the anisotropy ra- 
tios An = ilp,xx/^lp,*e; A VF = {vp^x} 1 / (vf,**) 1 f° r opti- 
mally doped LaCuO between a calculation in LAP W 51 ' 52 and 
the M31BM (model OP), respectively. fl p ,ij is given in units 
of THz and the Fermi velocity {vp^j} 1 ^ 2 in units of 10 7 cm/s. 




FIG. 6: Displacement patterns of the high-frequency oxygen 
bond-stretching modes (OBSM) of LaCuO. Left: apex-oxygen 
breathing mode (Of), middle: planar breathing mode (Og), 
right: half-breathing mode (Ai/2). 



strong softening of this mode across the insulator-metal 
transition was predicted prior to the experimental obser- 
vation, see Refsi 28 ' 39 i 42 . 

Optical modes like Of are of special importance for 
the cuprates because of their strong nonlocal coupling 
to the electrons due to the long-ranged Coulomb interac- 
tion poorly screened along the c-axis in these compounds, 
even in the metallic state. As a consequence nonlo- 
cal, nonadiabatic polar electron-phonon coupling effects 



become substantial. These high-frequency c-axis polar- 
ized optical phonons are basically unscreened. Quite re- 
cently such a long-ranged polar electron-phonon interac- 
tion (EPI) has been identified as very essential also for 
pairing in the cuprate superconductor o 53 i 54 . 

Displacements of the ions in the ionic layers of the 
cuprates, in particular those of the Of mode, bring about 
large changes of the potential felt by the electrons (holes) 
in the CuO plane which on her part are responsible for 
the superconductivity and the unusual normal state ef- 
fects. So, the c-axis polarized optical phonons become 
important for the renormalization of the electrons in the 
CuO plane, besides the renormalization due to strong 
correlation effects. Thus, both, the short-ranged and the 
long-ranged part of the Coulomb interaction is significant 
for the physics in the HTSCs. 

The manifestation of the coupling between electrons 
and polar c-axis phonons in the self-energy of the nodal 
quasiparticles in Bi2201 has been pointed out recently in 
Refs. [HEi. Note, that the long-ranged Coulomb cou- 
pling perpendicular to the CuO plane is very special for 
the HTSCs and would not be possible in a conventional 
three-dimensional metal or superconductor because of lo- 
cal screening by a high-density electron gas. As subse- 
quently discussed this long-ranged polar interaction does 
not only probe the in-plane electron properties but also 
leads to a strong coupling between electrons and c-axis 
polarized phonon modes. 

In the insulating parent compound Of is found at 
about 17 THz while in the optimally doped material 
the experimental value is estimated at about 11.5 THz 31 . 
However, the large linewidth of about 3-4 TH a 31 ' 32 and 
the presence of another Z-point mode at about 14 THz 
makes it difficult to pin down the exact position of Of 
in the metallic phase. 

In the following we will demonstrate within the 
M31BM that due to the nonlocal, nonadiabatic charge 
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response two coupled phonon-plasmon modes of Of -type 
arise at 9.09 THz and 13.86 THz, respectively, and that 
we can attribute to Of an adiabatic frequency of 12.96 
THz. Before discussing these results in more detail we 
present in retrospect a short explanation of the large soft- 
ening of Of during the insulator-metal transition, for de- 
tails see e.g. Ref. [28| and earlier references therein. 

In the insulating state Of can be shown to induce 
only an intralayer charge transfer between the copper 
and oxygen orbitals such that local charge neutrality of 
the cell is maintained under a perturbation due to Of. In 
this way the charge rearrangement leading to screening 
is considerably restricted and thus we find a high calcu- 
lated frequency of Of in the insulating state in agree- 
ment with experiment. Contrarily, in the metallic state 
no such restriction is present and Of generates CF's at 
Cu and O xy that have the same sign in the whole CuO 
layer. This ultimately leads to an (instantaneous) inter- 
layer charge transfer in adiabatic approximation which 
provides an effective screening mechanism for the long- 
ranged Coulomb interaction and produces the anoma- 
lous softening of Of during the insulator-metal transition 
in the calculations consistent with the measurements. 
Eventually, the instantaneous interlayer charge transfer 
in the adiabatic approximation is replaced by a collec- 
tive dynamic charge transfer within the phonon-plasmon 
scenario in a small region around the c-axis, including of 
course the Z-point. Such a scenario will be investigated 
next within the M31BM. 

In Fig. [7] we present our results of the dispersion of 
the coupled phonon-plasmon modes in the small region 
around the c-axis characterized by a nonadiabatic charge 
response. Specifically, we have calculated the coupled 
mode dispersion of the Aj modes from the Z-point along 
A' to the Z' = (e-^M, -tO point with e = 0.02, see 
Fig. [8l The mixed mode dispersion in Fig. [ZKa) is for 
the optimally doped case of LaCuO (model OP), Fig. 
EJb) and 7(c) show the results for the overdoped mate- 
rial, i.e. for model ODl and OD2, respectively. The 
strongly coupling Of modes are indicated as full dots. 
The broken line is the dispersion of the free-plasmon 
branch, calculated from Eq. (f2"3"| and the dotted line 
is the borderline for damping due to electron-hole decay, 
i.e. k ™g X z (E n (k) - E n (k + q)) for the band crossing the 
Fermi level. We observe a slight increase of the frequency 
of the free-plasmon branch upon doping which indicates 
a small reduction of the anisotropy. 

The phonon-plasmon scenario is exemplified for opti- 
mally doped LaCuO in Fig. [T^a) corresponding to the 
model OP for the proper polarization part. From in- 
spection of the figure we can learn two important facts. 
Firstly, the range of the region in q-space with a nona- 
diabatic charge response can be estimated and, secondly, 
the massive line broadening of the Of mode can be un- 
derstood. The full dots at the Z-point are the phonon- 
like and plasmon-like Of (na) mode at 9.09 THz and 
13.86 THz, respectively. We denote the lower mode as 
phonon-like because this mode connects to the adiabatic 
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FIG. 7: Calculated coupled phonon-plasmon dispersion of the 



A'i, (Of ) modes ( — ) from the Z-point (e = 0) along A' = 
(e^,0, 2zt) to £= 0.02 (Z') for the models OP (a), ODl (b) 
and OD2 (c), respectively. The coupling Of modes at Z are 

represented by •; : free plasmon branch; • ■ ■ borderline 

for damping due to electron-hole decay. 




FIG. 8: Schematic repre- 
sentation of the nonadi- 
abatic region in the (k x , 
k z ) plane with the direc- 
tions A' = (e^,0,2f), 
A' = £( e 2E,0,^) and 

A" = ( e ^,o,c%). <£ [0,1]. 
na: nonadiabatic region; ad: 
adiabatic region. 



Of (ad) phonon in the adiabatic region. Of (ad) has 
virtually the same frequency as the Of (ad) mode calcu- 
lated in adiabatic approximation, i.e. Of (ad) = 12.96 
THz. From the figure we can extract that the region of 
nonadiabatic charge response, characterized by the steep 
phonon-like branch is very small. It can be estimated at 
about e rj 0.01. 

The plasmon-like branch starts at the higher Of (na) 
mode, rapidly leaves the frequency range of the phonon 
spectrum within the nonadiabatic region and becomes 
more and more plasmon-like as it approaches the free- 
plasmon dispersion. 

The experimentally measured massive line-broadening 
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for Of of about 3-4 THa^^ in the optimally doped 
probe can be understood from the calculated phonon- 
plasmon scenario displayed in Fig. UJa). Experimentally, 
there is a limited wavevector resolution in the transverse 
direction perpendicular to the c-axis which is on aver- 
age e ~ 0.03 (Ref. HH). Thus, the relevant frequency 
range sampled by the neutron scattering experiment is 
over the steep branch and is between the two nonadia- 
batic Of (na) modes leading to an estimated broadening 
of about 3.9 THz. Because the region with a metallic, 
adiabatic charge response outweights by a factor of three 
the nonadiabatic region we can attribute to the broad 
Of mode an adiabatic frequency of about 12.9 THz. 

The calculated character of the coupled modes along 
A' can be obtained from Table HT1 where the displacement 
amplitudes for the phonon-like and plasmon-like mode 
is given in terms of relative displacements of the apex 
oxygen O z and the La ion along the c-axis. At the Z- 
point O z and La are vibrating in phase for the phonon- 
like mode and out of phase for the plasmon-like mode. 
In case of the phonon-like mode the amplitude of Of 
remains more or less constant while the amplitude for 
La decreases with increasing e and finally runs out of 
phase with O z in the adiabatic region. In the plasmon- 
like Of mode the O z and La ions always vibrate out 
of phase. Increasing e the amplitude of Of decreases 
while the magnitude of the amplitude of La increases. 
The different phase relation of the two modes is reflected 
in the corresponding induced charge redistributions as 
discussed below. 

The size e f» 0.01 of the nonadiabatic region predicted 
by our calculation within model OP is too small for an 
experimental observation of the dispersion of the coupled 
phonon-plasmon modes shown in Fig. [7j Only sampling 
the nonadiabatic region on average as mentioned above 
seems possible at present. 

Inspection of Fig. [7jb) and Fig. [7jc) for the overdoped 
states leads to about the same size of the nonadiabatic 
sector as in the optimally doped case and to the same adi- 
abatic frequency for Of (ad) and Of (ad), respectively. 
However, the broadening of Of is predicted by the cal- 
culation to decrease as compared to the optimally doped 
case. A decrease of the linewidth of Of with practically 
no change of the frequency has also been found quite re- 
cently for overdoped LaCuO in the experiments 57 . Thus, 
our prediction within the rigid-band approximation is 
conform with the measurements. From Fig. [3(b) and 
Fig. Etc) the width of Of is 3.23 THz for x w 0.22 
and 2.66 THz for x m 0.3, respectively. Last but not 
least we can read off from Fig. [7J a doping dependence 
of the phonon-like Of mode in which the mode energy 
decreases with reduced doping (OD2: 10.31 THz; OD1: 
9.74 THz; OP: 9.09 THz). 

Now we comment on the present interpretation of the 
measured dispersion of the steep Ai-branch displayed in 
Fig. [SJ We just have shown through our calculation 
of the phonon-plasmon dispersion along the A' direction 
that the nonadiabatic region around the c-axis is so small 



that it cannot be resolved by the INS experiments and 
only an average of the dispersion can be measured. So, we 
investigate in detail the coupled mode dispersion along 
the two directions A" = (e^,0, C^r) in Fig. Eland A = 

C { £ a>®>c) ^ n 10 nearby the c-axis in order to 

study the transition from the region of nonadiabatic to 
adiabatic charge response together with the change of the 
mode behavior during this transition. 

The frequency of the plasmon near V in Fig. [5] is al- 
ways very high and far outside the range of the phonon 
spectrum. Consequently the phonon dynamics is adia- 
batic in this q-space region. The phonon-like nonadia- 
batic Of mode (lower full dot) converges to its adia- 
batic value (open circle) at about e = 0.01, while the 
plasmon-like Of mode (upper full dot) rapidly leaves 
the range of the spectrum. Simultaneously the dynam- 
ics becomes adiabatic. Compare with the calculation of 
the phonon dispersion in the right most panel where the 
static (uj — 0) approximation for T1 KK > has been used. 
Likewise as in the investigations along the A' direction 
the nonadiabatic region can be characterized by e ~ 0.01. 
and the present resolution limit for INS does not allow to 
resolve the mixed-phonon plasmon dispersion. Only the 
result of averaging A" between T' and Z' can be detected. 
In such a procedure the small nonadiabatic part close to 
the c-axis will be outweighted in the measurement by a 
significant larger part where the dispersion is nearly adi- 
abatic. In particular we extract the development of the 
dispersion which leads to the steep Ai branch along A 
found in our calculations in the adiabatic limit. From 
this discussion it becomes clear how the measured data 
points in Fig. [2] should be interpreted, namely as an av- 
erage over the coupled mode dispersion with a dominant 
adiabatic contribution. 

Concerning the steep Ai branch the same conclusion 
can be drawn from the inspection of the calculated re- 
sults along the A' direction displayed in Fig. [TUJ Again 
the phonon-like Of mode matches the adiabatic result 
at about e w 0.01 and the dispersion with larger e val- 
ues is virtually adiabatic. At the T-point two coupling 
longitudinal A2 U modes (full dots) occur which we de- 
note as ferroelectric modes (FM) because of their intrin- 
sic displacement patterns where the oxygen anions are 
vibrating coherently against the cations in the lattice. 
As a consequence the electric dipole moments generated 
by the motion add constructively to a large value. So, we 
can expect a large oscillator strength and as seen in the 
experiments^ the FM dominates the infrared response 
for polarization along the c-axis not only in the insulat- 
ing state of LaCuO but also in the well doped metallic 
state. As already mentioned, such an optical activity in 
the metallic phase cannot be explained using the adi- 
abatic approximation for the calculation of the phonon 
dispersion as is usually done by applying static DFT for 
the metal, because there will be no LO-TO splitting (A2„ 
splitting in the present case) being a measure of the os- 
cillator strength in such a calculation. 

On the other hand, the observed infrared response in 
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Z(ad) Z(na) 


0.0001 


0.001 


0.002 


0.003 


0.004 0.007 


0.01 


0.02 


LO 


12.96 


9.09 


9.08 


9.27 


9.84 


10.59 


11.41 


12.61 


12.86 


12.94 


Oz(z) 


0.69 


0.66 


0.66 


0.67 


0.69 


0.70 


0.71 


0.69 


0.69 


0.69 


La(z) 


-0.15 


0.25 


0.25 


0.22 


0.14 


0.05 


-0.04 


-0.13 


-0.15 


-0.15 


LO 


12.96 


13.86 


13.86 


13.90 


14.02 


14.28 


14.81 


18.69 


24.38 




Z (Z) 


0.69 


0.68 


0.68 


0.68 


0.68 


0.67 


0.67 


0.61 


0.59 




La(z) 


-0.15 


-0.20 


-0.20 


-0.20 


-0.21 


-0.22 


-0.24 


-0.33 


-0.38 





Table II: Calculated displacement amplitudes of the phonon-like and plasmon-like O z modes in model OP along the A' = 
(e— , 0, 2zu direction between e — (zT-point) and e = 0.02. The frequency is given in units of THz. The left most data set is 
the result for Of in the adiabatic limit. The plasmon-like mode for e = 0.02 has not been calculated because of its frequency 
is far out of the relevant frequency range, see Fig. [TJa). 
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FIG. 9: Nonadiabatic coupled phonon-plasmon dispersion as calculated with model OP for the proper polarization part, repre- 
senting the optimally doped state of LaCuO from V' = (e^-, 0, 0) to Z' along the A" = (e— , 0, £— ) direction. In the different 
panels of the figure from left to right the following values" of e have been used: e = 0.003^0.004, 0.006, 0.01, 0.025, 0.05, 0.100. 
The right most panel shows the results of a calculation in adiabatic approximation for e = 0.100. Only the A" branches ( — ) 
coupling to the charge fluctuations are shown. : free-plasmon branch; • ■ • : borderline for damping. Dots at Z': 



o Of (ad), na: nonadiabatic; ad: adiabatic. 



Of (na) 



the metallic phase can be understood from our nona- 
diabatic results in Fig. [TO] Here we find a very large 
A 2u splitting for the longitudinal plasmon-like FM (up- 
per full dot at r) at 18.41 THz and the corresponding 
transverse FM (open circle at T) with a calculated fre- 
quency of 8.1 THz and 7.4 THz in the experiments^ 3 -, 
respectively. Simultaneously, there is a smaller splitting 
between the latter mode and the longitudinal phonon-like 
FM (lower full dot at T) at 5.80 THz. With increasing 
e, i.e. when approaching the region of adiabatic charge 
response, the LO-TO splitting is closed from below and 
the plasmon-like branch with the plasmon-like FM at T 
and the plasmon-like Of mode at Z' rapidly disappears 
out of the phonon spectrum. Moreover, the phonon dis- 
persion becomes adiabatic and the signature of the steep 
Ai branch appears in the A' direction. 

Summarizing, it should be emphasized that the pres- 
ence of the large nonlocal long-ranged polar nonadiabatic 



electron-phonon coupling which leads to the phonon- 
plasmon scenario is also reflected in the measured in- 
frared response. 

In Fig. [11] we show the calculated nonadiabatic 
phonon-plasmon dispersion of the Ai modes along the 
A ~ (0, 0, 1) direction. The full dots at T and Z represent 
the coupled phonon-like (lower dot) and plasmon-like A2 M 
FM (upper dot) and the phonon-like (lower dot) and 
plasmon-like (upper dot) Of mode, respectively. The 
calculated longitudinal plasmon-like A2 U FM at 18.41 
THz is not far away from the experimental result for the 
longitudinal A2„ FM of LaCuO in the insulating state 
obtained with INS at about 17.5 THz. Of course the 
large splitting between this mode and the correspond- 
ing TO mode (open dot) accounts for the dominance of 
the FM in the infrared response also in the insulating 
stat o 14 ' 43 . It should be remarked that in order to obtain 
the calculated dispersion curves the effect of the DF's is 
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FIG. 10: Nonadiabatic coupled phonon-plasmon dispersion as calculated with model OP for the optimally doped state of 
LaCuO from r to Z' = ?f) along the A' = C (e^.O, v) direction. In the different panels from left to right the 

following values for e have been used: e = 0.003, 0.004, 0.006?0.010 C , 0.025, 0.050, 0.100. The right most panel shows the result of 
a calculation within adiabatic approximation for e = 0.100. Only the A' branches ( — ) coupling to the charge fluctuations are 

shown. : free plasmon branch; ■ ■ ■ : borderline for damping. Dots at Z' : • Of (na), o Of (ad). Dots at V: • A^? (ferro, 

na); o A™ (ferro, ad), na: nonadiabatic; ad: adiabatic. 




FIG. 11: Calculated nona- 
diabatic coupled phonon- 
plasmon dispersion of the 
Ai modes of model OP for 
the optimally doped state of 
LaCuO along the A ~ (0, 0, 1) 
direction. The coupling 
modes at P (A^? (ferro, na)) 
and at Z (Of) are shown as 
black dots (•). The open dot 
(o) at P indicates A™ (ferro, 
ad). 

— : Ai modes; : free 

plasmon branch; ■ ■ ■ : border- 
line for damping; — ■ — : steep 
Ai branch in the adiabatic 
approximation. 
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important. The latter reduce the frequency of the plas- 
mon as compared with the calculation which takes only 
CF's into account. This effect is exemplified in Table 
IIIII for the free-plasmon and the plasmon- like modes at T 
and Z. 

From our preceding discussion of the coupled mode sce- 
nario along the A' and A" direction it can be concluded 
that the branch with the phonon-like FM at T and the 
phonon-like Of mode at Z' converges to the steep Ai 
branch of Fig. [IT] in the A direction when approaching 
the adiabatic limit. Thereby the LO-TO splitting at T is 
closed from below and the phonon-like Of mode reaches 
its adiabatic value. The plasmon-like branch connect- 



Table III: Frequency in THz of the free c-axis plasmon at the P 
and Z-point of the Brillouine zone (leftmost columns) as cal- 
culated with model OP for the optimally doped state. The no- 
tation CF means that only charge fluctuations are taken into 
account; CFD represents a calculation including additionally 
anisotropic dipole fluctuations. The rightmost columns gives 
the corresponding results for the plasmon-like modes. Com- 
pare with Fig. 1111 



ing the corresponding plasmon-like modes at T and Z' 
rapidly leaves the phonon-spectrum as can be seen in 
Fig. M 

A comment is appropriate regarding the relation be- 
tween electronic correlations in the CuO plane beyond 
LDA and the possible size of the nonadiabatic region 
which points to an interplay of nonadiabatic effects 
around the c-axis due to the long-ranged Coulomb in- 
teraction and strong short-ranged Coulomb repulsion re- 
sponsible for correlation. The latter leads in tendency to 
a reduction of the bandwidth and the Fermi velocity of 
the quasi-particle (QP) excitations in the (k Xl fc y )-plane 
as compared with DFT-LDA-like bandstructures^. Such 
a reduction should increase the size of the nonadiabatic 
region, because the electron dynamics in the CuO plane is 
slowed down and the free-plasmon frequency is accord- 
ingly decreased. In Ref. [59| we have simulated a nar- 
rowing of the electronic bandstructure within the 11BM 
by reducing empirically certain tight binding parameters 
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FIG. 12: Calculated nonadiabatic phonon dispersion of op- 
timally doped LaCuO in the main symmetry directions 
A ~ (1,0,0), E ~ (1,1,0) and A ~ (0,0,1) within 
the M31BM (model OP) for the proper polarization part. 
The various symbols representing the experimental results 
from Refs. I3lll32] indicate different irreducible representa- 
tion (ID's). Open symbols are for Lai.g Sro.i CUO4 and solid 
symbols for optimally doped LaCuO. The arrangement of the 
panels from left to right according to the different ID's is as 
follows: |Ai|A 2 (...,o), A 4 (— , □)|A 3 |E a |E 2 (-, o), E 4 (— ,D)| 
E 8 |Ai(— ,□), A 2 (- ■ ■ -,o)|A 3 |. 



and as a result a slight increase of the nonadiabatic re- 
gion around the c-axis has been found. Thus, it can be 
expected that characteristic properties of the QP in the 
CuO plane related to strong correlation effects take in- 
fluence on the dynamical charge response perpendicular 
to the CuO layer nearby the c-axis. 

Of course the position of the free-plasmon along the c- 
axis directly depends on the strength of the interlayer 
coupling which may differ in other cuprates from the 
situation found for LaCuO in the M31BM. As already 
mentioned, we have studied in the past parametrically 
the effect of a varying anisotropy in Refs. [li1fl"2l For a 
very weak interlayer coupling the resulting nonadiabatic 
dispersion in the metallic state along the A direction vir- 
tually cannot be distinguished from the dispersion of the 
(adiabatic) insulator and the optical phonon modes re- 
main completely unscreened. 

As far as the pairing channel via coupled phonon- 
plasmon modes is concerned it is interesting to remark 
that the total effective interaction between the electrons 
via phonon-plasmon exchange is attractive in the fre- 
quency range between the transverse A 2u FM and the 
longitudinal plasmon-like A 2u FM at the T point. Thus, 
the very large splitting of this mode, see Fig. [TTJland Fig. 
1 1 X I is favourable for pairing and the role played in this 
game by the incomplete dynamical screening of the (bare) 
polar Coulomb interaction in the nonadiabatic region be- 
comes evident. Note in this context that in recent work e 
it is found that the phonon-plasmon mechanism con- 
tributes constructively and significantly to the supercon- 
ductivity in the cuprates and that the long-range polar 
electron-phonon interaction is important for pairin g 53 ' . 
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FIG. 13: Same as in Fig. 1121 but calculated in the adiabatic 
approximation. 

In Figs. H31 we compare our calculations of the 
phonon dispersion of LaCuO along the main symmetry 
directions in the BZ within the M31BM as obtained for 
the nonadiabatic and adiabatic case of the proper polar- 
ization part, respectively. In Fig. [12] the small nona- 
diabatic region outside the c-axis where the adiabatic 
approximation grossly fails has been excluded because 
it cannot be resolved in this figure. Leaving aside the 
large differences of the calculated dispersion for the Ai 
modes as already discussed virtually no changes can be 
observed when passing from the nonadiabatic treatment 
to the adiabatic calculation. In particular this hold true 
for the modes propagating in the CuO plane. The reason 
is the fast electron dynamics in the CuO plane as com- 
pared to the dynamics perpendicular to the plane. From 
all the modes propagating in the CuO layer studied, the 
Og breathing mode, see Fig. O shows the largest de- 
viation, namely a stiffening of 0.284 THz (9.47 cm -1 ) 
for optimally doped LaCuO in the nonadiabatic treat- 
ment. We attribute this increase of the frequency to a dy- 
namical reduced nesting effect with wavevectors around 
X = ^(1, 1, 0). Such an interpretation is consistent with 
the calculated FS in Fig. 0Ja) and the maximum found 
at the A-point for the reduced electronic polarizability 
function n (q) from Eq. (|26|). The change of the FS 
topology upon doping from hole- like to electron-like and 
the related change of the nesting structures in Fig. 0] 
leads in the calculations to a reduction of the nonadi- 
abatic stiffening of Of (model OD1: 0.193 THz (6.44 
cm" 1 ), model OD2: 0.059 THz (1.97 cm^ 1 )). 

A possible effect of the nonadiabatic charge response 
on the Raman modes in LaCuO could be quite interest- 
ing because very recently giant nonadiabatic effects in 
layered metals such as graphite intercalated compounds, 
which are three-dimensional metals with a considerable 
anisotropy along one direction, have been predicted in 
Ref. Contrarily, in our calculations for LaCuO prac- 
tically no nonadiabatic effects have been found for the 
Raman modes. 

Finally, a comment on the calculated unstable mode in 
Figs. [TUnSlat the A-point should be made. This mode is 
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the so called tilt mode. Freezing of this distortion points 
correctly to the experimentally observed structural phase 
transition from the high-temperature tetragonal (HTT) 
to the low-temperature orthohombic (LTO) structure. 
In earlier calculations for the LTO structure we have 
found within the ab initio RIM two stable tilt modes 
with different frequencies in reasonable agreement with 
the experiment 12 . From these calculations it can be con- 
cluded that the structural phase transition is essentially 
driven by the strong long-ranged Coulomb interactions 
between the ions in the material which again emphasizes 
the important role of ionic binding for the cuprates. 



C. Phonon-plasmon induced charge response and 
selfconsistent changes of the potential 

In previous wor k 12 ' 61 ! 62 we have made explicit the 
charge response by calculating in adiabatic approxima- 
tion the phonon-induced charge density redistribution for 
the presumably generic phonon anomalies Ai/2 (half- 
breathing mode) and Og (planar-breathing mode), es- 
sentially polarized in the CuO plane and also for the Of 
mode 2 ^ polarized perpendicular to the plane, see Fig. 
[5] for the displacement patterns of these oxygen-bond- 
stretching modes (OBSM). We have shown in Sec. B that 
A i/2 and Og can be treated in the adiabatic limit, but 
this is of course not true for Of as a specific important 
mode from the nonadiabatic region of charge response. 

As suggested in Ref. [63| specific modes (in particular 
Ai/2) with a strong electron-phonon coupling can gener- 
ate an overscreening of the intersite Coulomb interaction 
due to a phonon-induced charge transfer between Cu and 
O xy . This strong coupling has been proposed in£2 to form 
a basis for the phonon mechanism of high-temperature 
superconductivity. 

The anomalous softening of the OBSM leading as in 
case of Ai/2 to a decreasing dispersion when starting at 
the r point, see Figs. 1121131 is a result of overscreening 
at shorter wave lengths of the changes of the Coulomb 
potential generated by the motion of the ions. This has 
been shown by our calculations to be due to nonlocally 
excited ionic CF's localized at the Cu and O xy sites gen- 
crating a dynamic charge ordering in form of localized 
stripes of alternating sign in the CuO layer. Such in- 
duced charge patterns also have been calculated quite 
recently for n-doped NdCuO 2 — . The importance for the 
phonon anomalies Ai/2 and Og of the repulsive on-site 
Coulomb interaction Ud as a relevant physical parameter 
of the cuprates deciding about the correlations in these 
materials has been predicted already in Ref. [H and in- 
vestigated in more detail in Refs. EgU Briefly, a large 
Ud suppresses double occupancy of the Cu3d orbital and 
thus Cu3d CF's which are crucial for the softening of Og 
and Ai/2. 

For the Ai/2 mode the charge stripes point along the 
CuO bonds and for Og along the diagonals of the CuO 
plane. Moreover, in Ref. [l2| we have pointed out by 



means of the OBSM, how antiferromagnetic correlations 
may favour the CF's in these modes and vice versa that 
nonlocal EPI is expected to generate antiferromagnetic 
spin-fluctuations (SF's) via the phonon-induced CF's be- 
tween the Cu ions. This interplay between charge- and 
spin-degrees of freedom should die out in the overdoped 
state because of the vanishing of the antiferromagnetic 
correlations. 

The induced charge redistribution for Of calculated in 
the adiabatic limit has been addressed shortly in Sec. B 
and will be compared explicitly in Fig. [T3]with the cor- 
responding results as obtained in the nonadiabatic ap- 
proach. So far we can say that the phonon anomalies 
Ai/2 and Og reflect the importance of the short-ranged 
part of the Coulomb interaction (Ud) while Of is specif- 
ically governed by the long-ranged part. 

Before proceeding with the investigation of the charge 
response in the optimally to overdoped metallic state of 
LaCuO we will include for completeness a short summary 
of our approach to phonon dynamics and charge response 
in the insulating and underdoped state of the cupra tes. A 

detailed discussion can be found in Refs. mumuM 

In this work also the strong doping dependence of the 
OBSM phonon anomalies has been studied. 

In order to describe these states of the cuprates we 
have proposed a microscopic modeling which takes into 
account the enhanced correlation of the Cu3d orbitals 
in these "phases" via a compressibilty-incompressibility 
transition for the Cu3d state in terms of rigorous sum 
rules for the proper polarization part II KK ' (q) in the 
longwavelength-limit. An interesting aspect of the par- 
tial incompressibility of the Cu3d states is that in case 
pairing correlations (pre-pairing) would exist above Tq 
the reduced local particle number fluctuations related to 
the orbital incompressibility of Cu3d due to the large 
on-site repulsion would correlate with enhanced phase 
fluctuations according to the particle number phase un- 
certainly relation and consequently a collective onset of 
phase coherence is suppressed. This in turn would stabi- 
lize a pre-paired state above Te- 
la Refs. 29.64 also the relation of the superconducting 
gap and the pseudogap as well as the possibility of Fermi 
surface transformations observed in the normal state of 
the cuprates is discussed within the model. Moreover, the 
modeling has been extended to n-doped cuprates and the 
electron-hole asymmetry is simulated via an orbital se- 
lective incompressibility-compressibility transition of the 
02p orbital in contrast to the p-doped case. 

We now return to the investigation of the phonon- 
induced nonadiabatic charge response in optimally doped 
LaCuO represented by the model OP. The nonlocal, non- 
rigid part of the charge response related to the nonlocal 
EPI effects as excited by a phonon with wavevector q 
and polarization a is given by 

5p n (r, qa) = £ a)p K (r - R*). (30) 

a.K 

The CF's S£* are obtained from Eq. (|19|) and the form- 
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FIG. 14: Contour plot in the CuO plane of the nonlocal part 
of the displacement-induced charge density redistribution for 
the Of mode in LaCuO calculated according to Eq. (|30[) m 
adiabatic approximation for the M31BM (model OP) for the 
proper polarization part (a), (b) same as in (a) but for the 
phonon-like Of mode, (c) same as (a) but for the plasmon- 
like Of mode. The units are in 10 -4 e 2 /a%. The phase of 
Of is as shown in Fig. [6] Full lines ( — ) mean that electrons 
are accumulated in the corresponding region of space and the 
broken-dot lines (—•••—) indicate regions where the electrons 
are pushed away. 



<1,0,0> 

(c) 



factors p K from Eq. 

In the contour plots of Fig. Q3] the displacement in- 
duced rearrangement of the charge density according to 
Eq. ([50)) is shown for the Of mode. The calculation 
has been performed for the model OP. Fig. [P4Ta) dis- 
plays the result for the adiabatic approximation and Fig. 
H~4T b) and Fig. [l"4tc) pictures the nonadiabatic result 
for the phonon-like and plasmon-like mode Of mode, re- 
spectively. The phase of the vibration of the ions is as 
in Fig 6, i.e. the apex oxygens move away from the CuO 
layer. The full lines (Sp n > 0) indicate that electrons 
are accumulated in the associated region of space. The 
broken lines (8p n < 0) mean that electrons are depleted 
in this region. 

As already mentioned in Sec. B and displayed explic- 
itly in Fig. 14(a) in the adiabatic approximation Of 
(12.96 THz) generates CF's at the Cu and O^y ions of 
the same sign in the whole CuO layer. In the adjacent 
layers the induced CF's are the same in magnitude but 



of opposite sign, compare with Fig. [PoT a). Thus, Of 
activates an interlayer charge transfer which is instan- 
taneous in adiabatic approximation. Such a nonlocally 
excited screening process is very effective for screening 
the changes of the Coulomb interaction generated by the 
displacement of the apex oxygen ions and explains ac- 
cordingly the large softening of Of during the insulator- 
metal transition. 

Fig. [PfTb) records our calculated results for the nona- 
diabatically induced CF's of the phonon-like Of mode in 
the CuO plane at 9.09 THz. In the nonadiabatic regime 
the electrons cannot follow instantaneously the move- 
ment of the ions. The now dynamically excited CF's 
at Cu and O xy have the same sign in the CuO layer. The 
sign is as in the adiabatic case because the frequency of 
the free plasmon at Z is (slightly) higher (9.5 THz, see 
Table IIII[) than the phonon-like Of mode. The magni- 
tude of the CF's is considerably larger in the nonadia- 
batic case, see also Table HVl 
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FIG. 15: Contour plot perpendicular to the CuO plane of 
the nonlocal part of the displacement-induced charge density 
redistribution for the Of mode in LaCuO calculated from Eq. 
(I3U|) in adiabatic approximation for the M31BM (model OP) 
(a), (b) same as (a) but for the phonon-like Of mode, (c) 
same as (a) but for the plasmon-like Of mode. The units 
and the meaning of the line types as in Fig. 1141 The arrows 
denote the displacements of the apex oxygen ions. 



The interlayer CT is now of dynamic nature because 
of the dynamical screening of the Coulomb interaction 
and consequently a coupled collective phonon-like exci- 
tation propagating along the c-axis is created leading to 
the correlated charge rearrangement between the layers 
due to the coherent dynamics as shown in the snapshot of 
Fig. [ToT b). Comparing with the adiabatic limit besides 
the strongly enhanced charge response, which explains 
the lower frequency, the sign of the CF's at the La ion 
has changed as compared with the adiabatic calculation. 
This correlates well with the change of sign of the dis- 
placement amplitude for La, compare with Table [TTJ On 
the other hand, the sign of the CF's at the La for the 
plasmon-like mode is the same as in the adiabatic case 
in agreement with the same sign of the displacement of 
La in both calculations. 

The coupled plasmon-like Of mode at 13.86 THz has 
a higher frequency than the free plasmon and appar- 
ently the charge response runs out of phase, see Fig. 
I14f c). This results in an antiscreening consistent with 
the higher mode frequency. The sign of the CF's at Cu 
and O xy is not only opposite to that in the adiabatic and 
phonon-like case but also smaller in magnitude (Table 
ITV|) . Fig. [Mtc) together with Fig. [TBTc) give a snapshot 
of the collective plasmon-like interplane charge excita- 
tions being out of phase with the lattice and also strongly 
reduced in magnitude as compared with the phonon-like 
excitation displayed in Fig. [I4jb) and Fig. [TM b) , respec- 
tively. 




FIG. 16: Calculated displacement-induced charge-fluctua- 
tions 5f re (q, a, ui) (Eq. (|19|l ) (a) and corresponding orbital 
averaged changes of the crystal potential 5V K {<\, cr, u) (Eq. 



22)) 1 (b) for the phonon-like Of mode along the A' direc- 



tion (e^,0, f),e£ (0.00,0.02). Cu3d (— 
orbital degrees of freedom are considered, 
have been performed for model OP. 



I and O xy 2p (• ■ ■ ) 
The calculations 



Concerning c-axis transport in a highly anisotropic 
metal it is interesting to note that in a simplified model 
treatment the possibility has been pointed out that a 
strong coupling between electrons and a bosonic mode 
propagating and polarized in the c-direction, like Of, 
can lead to a metallic to nonmetallic crossover in the c- 
axis resistivity as the temperature is increased with no 
corresponding feature in the xy plane properties 65 . 

In order to get an impression for the displacement in- 
duced changes of the CF's <5Ck(cu cr) according to Eq. (TT5)) 
and the corresponding induced orbital averaged changes 
of the potential felt by the electrons 8V*(q, a) from 
Eq. (|22[) , when passing from the region with a nonadia- 
batic dynamic charge response to the region with practi- 
cally static adiabatic charge response, we display in Figs. 
ini[a,b) our calculated results for the phonon-like Of 
mode along the A' direction, compare with Fig. EJa). 
The calculations have been carried out for the optimally 
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doped case using model OP and for the most important 
Cu3d and O xy 2j) orbitals, respectively. Similar as for the 
phonon frequencies in Fig. [7ja) S( K and 8V K converge be- 
yond e w 0.01 to their values in the adiabatic limit. Tak- 
ing for example the calculated numbers at e = 0.00 and 
e = 0.02, respectively, we recognize that the nonlocally 
excited CF's and in particular the nonlocally induced or- 
bital averaged changes of the potential being a measure 
of the long-ranged polar coupling of the electrons and 
phonons are very strongly enhanced in the nonadiabatic 
region of phase space around the c-axis. 

Additionally, in Table IIVI we summarize and compare 
the calculated orbital selective values of <5£k an d o~V K for 
Cu3d and O a , l) 2p of the phonon anomalies Ai/2, Og and 
Of, respectively. We extract already strong coupling of 
the adiabatic modes Ai/2, Og and of course the strong 
enhancement for the phonon-likc Of mode. This again 
points to the fact that besides the short-ranged part of 
the Coulomb interaction Ud being important for Ai/2 
and Og, the strong long-ranged, nonlocal polar Coulomb 
interaction must play an important role in the normal 
and superconducting state of the doped cuprates. 

In context with the very strong coupling of the doping 
dependent phonon-like Of mode which is at 9.09 THz 
(37.6 meV) in model OP for the optimally doped state 
a special remark should be made. This mode involves 
momentum transfer between antinodal regions of the FS 
(similar as Og or the resonant magnetic mode), how- 
ever, with a drastically enhanced coupling strength as 
compared to Og , see Table [IVJ Consequently, this mode 
should contribute to the electron self-energy via electron- 
phonon coupling in the antinodal region and possibly also 
to the antinodal pseudogap in the normal state. In the 
superconducting state even an increased contribution to 
the self-energy via Of can be expected because of the 
large density of states enhancement in these regions due 
to the opening of a gap with d-wave symmetry^. More- 
over, below Tc we will have a suppression of the low scat- 
tering rate of the electrons due to the opening of a gap 
and additionally below the finite energy of the phonon- 
like mode. 

This might be helpful to understand the dramatic 
change in lineshape of the antinodal spectra seen in 
ARPES, where a sharp quasi-particle peak develops at 
the lowest binding energies followed by a dip and a 
broader hump. In the current literature these features 
are explained by a finite energy collective bosonic mode 
which mostly is identified as the magnetic resonance 
mod o 67 ' 68 . On the other hand, the enhanced coherence 
of the quasi-particles at lower temperature in the super- 
conducting state is simultaneously favourable for the ex- 
istence of a phonon-plasmon scenario with the plasmon- 
like mode as another possibility for the collective bosonic 
mode. 

From our calculations damping of the plasmon due 
to electron-hole decay is less important because only in- 
traband excitations contribute in case of the simple BS 
near the Fermi level. So, the phonon-like mode should 



be an important low-energy bosonic excitation mediating 
an attractive interaction. Comparing with the isostruc- 
tural low-temperature superconductor SrRuO these re- 
lations can be expected to change considerably because 
of the complex BS around the Fermi energy where three 
FS sheets exis t 69 i 70 ' 71 . In such a situation a strong q- 
dependent plasmon damping by low-lying interband tran- 
sitions can be expected in the collisionless regime and the 
c-axis plasmon might cease to be a well defined excita- 
tion. Note, that in simple metals is has been shown that 
interband transitions are the dominant processes in the 
damping of plasmons 72 . Additional damping due to in- 
teractions between electrons and with other degrees of 
freedom enhances the damping still further and broad- 
ens the energy levels. This leads to a narrowing of the 
interband gaps and to the fact that the decay via inter- 
band transitions would even be promoted. Altogether, 
it seems likely that in contrast to LaCuO in the low- 
temperature superconductor SrRuO the c-axis phonon- 
plasmon modes are not well defined. On the other hand, 
the existence of these modes in the cuprates could point 
to their relevance for high-temperature superconductiv- 
ity. 

In case of cuprate superconductors with more than one 
CuO layer a qualitatively new possibility may arise as 
candidate for the finite energy collective bosonic mode. 
In these compounds additionally to the low-lying plas- 
mon (intraband plasmon) of the single layer compounds 
discussed so far which results from a weak but non van- 
ishing k z dispersion of the QP's an interband plasmon 
with a low energy can exist. The latter emerges at van- 
ishing k z dispersion from a finite but small intra-layer 
coupling, e.g. between the two CuO layers in a double 
layer system. Such a possibility has been investigated 
qualitatively for YBaCuO within a simplified model ap- 
proach, see Refs. JT2JI73I analogously to the model for 
LaCuO in Refs. IIOGJ. 

In case of YBaCuO a 22BM is used and then the de- 
generacy of the BS is lifted by introducing parametri- 
cally an intra-bilayer coupling. The latter is varied to 
simulate a corresponding bilayer splitting. According 
to these calculations a sufficiently weak bilayer coupling 
generates a nearly dispersionless low-lying phonon-likc 
bond bending branch polarized along the c-axis with a 
phonon-like (O23 TT) mode at the Z-point and a lon- 
gitudinal phonon-like FM (B{ u (0 2 3 TT)) at the T-point. 
The notation (O23 TT) means that the O2 and O3 oxy- 
gen ions in the CuO layer are vibrating in phase. With 
increasing bilayer-coupling the frequency of the phonon- 
like Ag (O23 TT) varies between around 35 meV at weak 
coupling and 46 meV in the adiabatic limit. 

The calculated phonon-induced charge rearrangement 
of A z g (0 23 TT) m YBaCuO is similar to that of Of in 
LaCuO as far as a double layer of YBaCuO is concerned. 
The CF's have the same sign in the CuO plane but are of 
opposite sign in the two CuO planes in the double layer. 
This is energetically favourable because of the attractive 
Coulomb interaction at short distance between the CF 
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Of (ad) 


Of(na,ph) 


Of (na.pl) 


Of (ad) 


Ai/2(ad) 


Of (na) 


Ai/2(na) 




Cu3d 02p 


Cu3d 02p 


Cu3d 02p 


Cu3d 02p 


Cu3d 02p 


Cu3d 02p 


Cu3d 02p 


sv K 


-6.59 -2.63 
-15.77 58.67 


-34.12 -18.13 
-955.98 -515.98 


1.32 2.13 
241.34 207.04 


-16.45 
-93.92 


-14.76 -1.44 
-73.06 55.14 


-16.95 
-87.20 


-14.76 -1.44 
-73.13 55.08 



Table IV: Magnitudes of the charge fluctuations 5£« (Eq. (|19p ) and orbital averaged changes of the selfconsistent changes of 
the crystal potential 8V K (Eq. (|22[0 of the OBSM as calculated with model OP for the Cu3d and the 02p orbital in the CuO 
plane. SCk is given in units of 10 -3 particles and 8V K in units of meV. na: nonadiabatic; ad: adiabatic; ph: phonon-like; pi: 
plasmon-like. The phase of the displacements for the OBSM is as shown in Fig. H3 SCk < means an accumulation of electrons 
in the corresponding orbital. SV K > indicates that the region around the ion is attractive for electrons. 



redistribution of the two layers. So, this CF pattern also 
explains the strong renormalization of the frequency of 
this mode as compared with a calculation where the CF's 
and the related coupling within de bilayer is ignored (61 
meV) and DF's only are considered^. 

Concerning the infrared response the weak bilayer- 
coupling only leads to a partial screening of the long- 
ranged Coulomb interaction and thus to a corresponding 
fractional reduction of the LO-TO splittings, i.e. the lat- 
ter are still present consistent with the measured optical 
c-axis spectra in YBaCuO, see e.g. Refs. IzB How- 
ever, the splittings are gradually closed from below by 
admitting additionally to the intra-bilayer coupling an 
inter-bilayer coupling in the model leading to a non van- 
ishing k z dispersion and, finally, to the complete metallic 
screening of the long-ranged Coulomb interaction in the 
adiabatic limit. 

The Ai/2 anomaly involves momentum transfer be- 
tween nodal regions and could be important to under- 
stand the corresponding self-energy corrections, in par- 
ticular, the nodal kink observed in ARPES experiments. 
Indeed, in a recent bulk-sensitive low photon energy 
ARPES studySi it has been demonstrated that electron- 
phonon coupling is responsible for the nodal kink and the 
"half-breathing" mode Ai/2 is the relevant mode. 

Another general aspect becomes evident from our cal- 
culations namely that the specific features of the solid- 
state chemistry of the cuprates are essential. The dis- 
tinctive ionic character of these compounds underlines 
the relevance of the long-ranged polar Coulomb inter- 
action also in the doped metallic state. On the other 
hand, the Mott insulating character of the undoped par- 
ent compounds being antiferromagnetic insulators em- 
phasizes the importance of the short-ranged part of the 
interaction, notably by the repulsive on-site interaction 
Ud- As discussed in this work the impact of both parts of 
the interaction appears when studying the generic OBSM 
phonon anomalies. 

In the overwhelming number of attempts to describe 
the physics of the cuprates the short-ranged part Ud 
is considered exclusively and simplified models with 
strongly reduced Hamiltonians are used. Such models in 
the first place may be useful for the discussion of theoret- 
ical ideas in general rather than for realistic calculations 
where the material specifica must be included. Thus, 
at present it is not at all clear if the two-dimensional 



one-band Hubbard model which is most commonly ap- 
plied for the cuprates contains sufficient information to 
describe the real materials. In this context a recent com- 
bined DFT and dynamical cluster Monte Carlo calcula- 
tion for the three-band Hubbard model is instructive and 
may serve as a guideline^ in the discussion of a minimal 
model for the cuprates. This study demonstrates that 
the corresponding phase diagram of the model and in 
particular a possible occurrence of superconductivity is 
quite sensitive to the choice of hopping parameters and 
the down-folding procedure used. In a recent study££ 
pairing correlations on doped Hubbard models are rean- 
alyzed using a sign-problem-free Gaussian-Basis Monte 
Carlo method and it is concluded that the simple Hub- 
bard model does not account for high-temperature su- 
perconductivity. 



The strong short-ranged repulsive Ud by itself does not 
help with the pair-binding, however, favours d-wave sym- 
metry and the usual argument then is that as an implica- 
tion of a large Ud an antiferromagnetic exchange coupling 
J arises which leads to an attractive interaction between 
electrons of opposite spins on neighbouring sites. Ud is 
an unretarded particle-particle interaction with no low- 
frequency dynamics and also the interaction related to J 
is unretarded. 



On the other hand, several retarded attractive interac- 
tions mediated by virtual bosonic excitations with low- 
frequency dynamics are discussed in context with pair- 
ing in the cuprates. Most common are fluctuations of 
spin and/or charge degrees of freedom of the electron liq- 
uid and (adiabatic) phonons. From our findings in this 
paper low-frequency coupled phonon-plasmon modes re- 
sulting from the strong polar long-ranged interaction in 
the nonadiabtic region should be added to the list. It 
still remains an open question which contribution and 
what kind of cooperation between the different players is 
most important for pairing and superconductivity in the 
cuprates. It seems that the high as well as the low energy 
scale is essential for a synergetic interplay of charge-spin- 
and lattice degrees of freedom which underlies the physics 
in the cuprates in the normal as well as in the supercon- 
ducting state. 
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IV. SUMMARY AND CONCLUSIONS 

We have developed a realistic description of the elec- 
tronic bandstructure of LaCuO which is well suited to 
calculate the wavevector and frequency dependent proper 
polarization part of the DRF for the optimally to over- 
doped state. The latter is used to calculate the adiabatic 
and nonadiabatic charge response and the coupled mode 
dynamics in LaCuO. 

The large anisotropy along the c-axis of the electronic 
structure of the cuprates is considerably underestimated 
in DFT-LDA calculations. So we have modified a LDA- 
based 31BM to account for the much weaker interlayer 
coupling in the real material. We have optimized the 
interlayer coupling by reducing some relevant hopping 
parameters perpendicular to the CuO layer in such way 
that significant features of the Ai phonons polarized par- 
allel to this axis are well described. The Ai modes have 
been chosen because they are most sensitive in respect to 
the charge response orthogonal to the CuO plane. The 
new sufficiently anisotropic model (M31BM) leads to a 
more extended saddle point region, a reduction of the 
width of the bandstructure, flatter bands in c-direction 
and a significant narrowing and increase of the important 
van Hove peak. Simultaneously, we obtain an associated 
amplification on the average of the proper polarization 
part. These results demonstrate the importance of a cor- 
rect representation of the interlayer couplings being rele- 
vant for both the electronic properties in the CuO layers 
and along the c-axis. 

Within the M31BM the measured Fermi surfaces for 
optimally and overdoped LaCuO are well described. 
Clearly visible in the experiments and the calculations is 
the change of the FS topology and of the nesting struc- 
tures upon doping. The latter are shown to be reflected 
in the (noninteracting) susceptibility. Moreover, relevant 
Fermi surface parameters for transport like the Drude 
plasma energy tensor and the Fermi velocity tensor have 
been calculated and compared with standard LAPW cal- 
culations. An enhancement of about a factor of 5 for the 
anisotropy ratio is found for both types of parameters in 
the M31BM. 

In our calculation of the phonon dispersion of the c-axis 
polarized Ai modes dramatic changes occur between the 
adiabatic limit based on a static DRF and the nonadia- 
batic calculation founded on a dynamic DRF. The former 
leads to static and the latter to dynamic screening of the 
Coulomb interaction. On the other hand, virtually no 
nonadiabatic effects have been detected for the phonons 
propagating in the main symmetry direction A and X in 
the CuO plane. This is due to the fast dynamics of the 
electrons in the CuO layer as compared with the slow 
dynamics perpendicular to the plane. Only the planar 
oxygen breathing mode Og experiences a minor nona- 
diabatic correction worth mentioning which we attribute 
to dynamical reduced nesting. 

The adiabatic approximation as a basis paradigma to 
investigate lattice vibrations fails severely in a small sec- 



tor of q-space around the A direction. The large nonadi- 
abatic effect found for the Ai modes results from a low- 
lying c-axis plasmon predicted within the M31BM. This 
collective mode couples via the strong long-ranged polar 
Coulomb interaction to the optical phonons of allowed 
symmetry. In particular we have shown that the large 
softening of the Of mode during the insulator-metal 
transition and its giant linewidth found experimentally 
can be well understood in the phonon-plasmon scenario. 

At the r point we identify ferroelectric-like Ai u modes 
with a large LO-TO splitting in form of plasmon-like 
and phonon-like excitations. As seen in the experiments 
the ferroelectric-like Ai u mode dominates the infrared re- 
sponse for polarization along the c-axis not only in the 
insulating state of LaCuO but also in the well doped 
metallic state. Such an optical activity also found in 
our nonadiabatic calculations for the metallic phase can- 
not be explained using the adiabatic approximation with 
static screening for the calculation of the phonon disper- 
sion. This is what has been routinely done in the litera- 
ture in first principles calculations within static DFT. So, 
the latter are not consistent with experimental evidence 
concerning the c-axis response. 

From our systematic investigation of the nonadiabatic 
charge response and phonon dispersion we can identify a 
small region in q-space around the A direction where the 
adiabatic approximation breaks down. In this sector of 
q-space the strong nonlocal, nonadiabatic polar electron- 
phonon coupling leads to a mixed phonon-plasmon dis- 
persion. In order to experimentally resolve the calcu- 
lated dispersion a very high resolution perpendicular to 
the c-axis would be needed. Presently, inelastic neutron 
scattering experiments with such a high resolution arc 
out of range and the actual measurements only perform 
some average orthogonal to the c-axis where the adia- 
batic charge response dominates by far because of the 
smallness of the nonadiabatic zone. 

We have calculated the induced charge response of the 
apex oxygen stretching mode Of in the adiabatic limit 
and the nonadiabatic phonon-plasmon regime. Addition- 
ally we have computed the associated orbital averaged 
changes of the potential felt by the electrons when pass- 
ing from the nonadiabatic phonon-plasmon region to the 
adiabatic regime. In the adiabatic case we have an in- 
stantaneous interlayer charge transfer while in the nona- 
diabatic situation the interlayer charge transfer is dy- 
namic in nature and phonon-like and plasmon-like exci- 
tations are excited. Comparing with the adiabatic result 
we find for the phonon-like Of mode a strongly enhanced 
in phase charge response which may be denoted as over- 
screening with respect to the adiabatic case and accord- 
ingly a lower frequency appears while for the plasmon-like 
mode with a higher frequency the charge response is sig- 
nificantly weaker and runs out of phase, i.e. we observe 
an antiscreening effect. 

The related changes of the orbital averaged potential 
have been calculated for the coupled phonon-plasmon 
modes. In particular for the phonon-like mode we have 
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found a very strong enhancment in the nonadiabatic sec- 
tor of phase space as compared with the adiabatic result. 
The possible importance of the phonon-like Of mode 
for the electron self-energy has been pointed out. This 
demonstrates again the importance of the long-ranged 
polar coupling of the electrons and the c-axis phonons. 
As far as the low energy scale and the bosonic glue to 
keep the electron pairs in the cuprates together is con- 



cerned our calculations make a strong case that phonon- 
plasmon modes from the nonadiabatic zone around the 
c-axis have to be added to the list of possible players. Fi- 
nally, it has been speculated within the phonon-plasmon 
scenarium why SrRuO could be a low-temperature su- 
perconductor in contrast to the high-temperature super- 
conductor LaCuO. 
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